Mass-Energy and Momentum Extraction by Gravitational Wave Emission in the 
Merger of Two Colliding Black Holes: The Non-Head-On Case 
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We examine numerically the post-merger regime of two Schwarzschild black holes in non head-on 
collision. Our treatment is made in the realm of non-axisymmetric Robinson- Trautman spacetimes 
which are appropriate for the description of the system. Characteristic initial data for the system 
are constructed and the Robinson- Trautman equation is integrated for these data using a numerical 
code based on the Galerkin spectral method which is accurate and sufficiently stable to reach the 
final configuration of the remnant black hole, when the gravitational wave emission ceases. The 
initial data contains three independent parameters, the ratio mass a of the individual colliding 
black holes, the boost parameter (that characterizes the initial pre-merger infalling velocity of the 
two black holes) and the incidence angle of collision < po < 90°. The remnant black hole is 
characterized by its final boost parameter, the final rest mass and scattering angle. The motion 
of the remnant black hole is restricted to the plane determined by the directions of the two initial 
colliding black holes, characterizing a planar collision. The energy-momentum fluxes carried out by 
gravitational waves are confined to this plane. We evaluate the efficiency of mass-energy extraction, 
the total energy and momentum carried out by gravitational waves and the momentum distribution 
of the remnant black hole for a large domain of initial data parameters. Our analysis is based on 
the Bondi-Sachs four momentum conservation laws. The process of mass-energy extraction is shown 
to be less efficient as the initial data departs from the head-on configuration. Head-on collisions 
(po = 0) and orthogonal collisions (po = 90°) constitute, respectively, upper and lower bounds to the 
power emission and to the efficiency of mass-energy extraction. On the contrary, head-on collisions 
and orthogonal collisions constitute, respectively, lower and upper bounds for the momentum of 
the remnant. The momentum extraction and the pattern of the momentum fluxes, as a function 
of the incidence angle, are examined. The momentum extraction characterizes a regime of strong 
deceleration of the system. The angular pattern of gravitational wave signals are also examined for 
small mass ratios a and early times u. We show that, in the plane of collision {x,z), the pattern 
is typically bremsstrahlung, corresponding to a strong deceleration regime at early times, with two 
dominant lobes in the forward direction of motion. Gravitational waves are also emitted outside 
the plane of collision, the angular distribution of which is symmetric with respect to the y-axis and 
therefore with a zero net momentum, consistent with the planar nature of the collision. Finally the 
relation between the incidence angle and the scattering angle closely approximates a relation for the 
inelastic collision of classical particles in Newtonian dynamics. 
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previous recent papers [H, Q where we examined the case 
of head-on collisions. The collision and merger of two 
black holes are considered to be important astrophysi- 
cal sources of strong gravitational emission (cf.[J| and 
references therein) and the relevance of the processes of 
generation and emission of gravitational waves in these 
configurations lies in the fact that the associated wave 
patterns will be crucial for the present efforts towards a 
direct detection of gravitational waves. In spite of the 
enormous progress achieved until now using approxima- 
tion methods and numerical techniques, the information 
on wave form patterns and radiative transfer processes in 
the dynamics of gravitational wave emission is far from 
being complete 5]. 

The object of this paper is the numerical treatment 
of the gravitational wave production and the associated 
processes of radiative transfer in the non-head-on colli- 

sion of two boosted Schwarzschild black holes, modeled 

in the context of Robinson- Trautman (RT) spacetimes. 

* Electronic address: ''r afael.aranha®physics.gatech.edu;ivanoQcbpf.br;tofaitlBaaiiffeeaa3tip initial data for the RT dynamics are COn- 



I. INTRODUCTION 



It is by now theoretically well established that, in 
the nonlinear regime of General Relativity, gravitational 
waves extract mass, momentum and angular momen- 
tum of the source and that the radiative transfer in- 
volved in these processes may turn out to be funda- 
mental for the astrophysics of the collapse of stars, the 
formation of black holes and the collision and merger 
of two or more black holes. In the present paper we 
will examine the radiative transfer processes by gravi- 
tational wave emission in the post-merger phase of two 
colliding black holes, with focus on the case of non- 
head-on collisions, described in the realm of Robinson- 
Trautman (RT) spacetimes This work extends two 
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structed that represent instantaneously two black holes 
in non-head-on collision. RT spacetimes present a global 
apparent horizon Q so that the dynamics corresponds to 
a regime where the merger has already set in. During 
the merger part of the rest mass and kinetic energy of 
the two initial individual black holes are radiated away 
by the gravitational waves, and part will be absorbed to 
constitute the mass-energy of the remnant. The linear 
momentum of the initial system is also carried out by 
the gravitational waves emitted. We also examine the 
two regimes of momentum extraction in the post-merger 
phase by contemplating the time behavior of the momen- 
tum flux of the gravitational waves emitted. We use the 
momentum-energy conservation law in the Bondi-Sachs 
formulation of gravitational wave emission by bounded 
sources 0, Also the analysis of the energy flux car- 
ried out by gravitational waves will allow us to charac- 
terize two distinct regimes of bremsstrahlung emission for 
distinct domains of the initial data parameters. Several 
characteristics of a non-head-on collision are drastically 
distinct from a head-on collision, mainly for large values 
of the initial mass ratio parameter. The general outcome 
will be a remnant boosted black hole with a rest mass 
larger than the sum of the rest masses of the two indi- 
vidual initial black holes. 

As compared to previous estimates of the literature, 
our approach differs basically in that we have adopted 
the characteristic surface initial data formalism, which 
has several advantages for the description of gravitational 
radiation and the construction of algorithms [9|j. In this 
direction, an accurate code based on the Galerkin method 
was constructed to integrate the RT non-axisymmetric 
field equations. The code is accurate and highly sta- 
ble for long time runs in the nonlinear regime so that 
we are able to reach numerically the final configuration 
of the system, when the gravitational emission ceases. 
This will allow us to describe processes of mass-energy, 
and momentum extraction due to gravitational radiation 
emission, by using physically meaningful quantities con- 
nected to initial and final configurations of the source. 

We organize the paper as follows. In Sec. II we re- 
view some basic aspects of the general non-axisymmetric 
Robinson- Trautman spacetimes necessary for our anal- 
ysis of gravitational radiation emission in the case of a 
non-head-on collision. In Section III we construct char- 
acteristic initial data for the non-axisymmetric RT dy- 
namics corresponding to two boosted Schwarzchild black 
holes in a non-head-on collision. In Section IV we present 
a basic description of the numerical codes used to evolve 
the characteristic initial data via the RT dynamics. The 
codes are based on Galerkin spectral methods and the 
dynamical evolution of the initial data is performed and 
discussed, including its accuracy and stability. In Sec- 
tion V we discuss the planar nature of the collision and 
show how this property allows us to save a lot of com- 
putational effort. In Section VI the Bondi-Sachs four 
momentum for the non-axisymmetric RT spacetimes are 
introduced together with its conservation laws that will 



be fundamental in the treatment of the energy and lin- 
ear momentum extracted from the source by the gravita- 
tional waves emitted. Sections VII to X contain the main 
results of the paper, as the numerical evaluation of the 
efficiency of mass-energy transfer by gravitational wave 
emission, the conservation of energy in the processes via 
the Bondi mass formula and linear momentum extrac- 
tion by gravitational waves. Distinct patterns of gravi- 
tational wave emission are discussed related to the range 
of the physical parameters of the initial data. In Section 
XI we summarize our results and discuss their relevance 
and limitations as compared to previous results in the 
literature. Throughout the paper we use units such that 
87rG = c= 1. 



II. DYNAMICS OF ROBINSON-TRAUTMAN 
SPACETIMES. WAVE ZONE CURVATURE 

RT spacetimes [l| are asymptotically fiat solutions of 
Einstein's vacuum equations that describe the exterior 
gravitational field of a bounded system radiating gravi- 
tational waves. The RT metric is expressed as 



ds =Q {u,r,6,(t>)du +2dudr~r K {u,0,(, 

X {de^ + sin edif^ 



(1) 



where r is an affine parameter defined along the shear- 
free null geodesies determined by the vector field d/dr. 
Einstein equations imply that, in a suitable coordinate 
system, 
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^X{u,eA)- — + 2r^^., (2) 
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where uiq is a positive constant, and A(u, 0, 0) is the 
Gaussian curvature of the surfaces (u^const, r=const) 
defined by 
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The remaining Einstein equations yield 
'(Aesin6')e A^^ 
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= 0. 



(4) 



In the above, subscripts u, 9 and </> denote derivatives 
with respect to u, 9 and </>, respectively. Eq. (|H), de- 
noted RT equation, governs the dynamics of the grav- 
itational field which is totally contained in the metric 
function K{u,9,(t)). Chrusciel and Singleton pLO,] estab- 
lished that RT spacetimes exist globally for all positive u 
and converge asymptotically to the Schwarzschild metric 
as u — oo - this global time extension being realized for 
arbitrary smooth initial data. 

An important feature of RT spacetimes, that estab- 
lishes its radiative character, arises from the expression 
of its curvature tensor that in the semi-null tetrad basis 



= (qV2) du + dr 



e" = du, 



= rKde, Q-^ = ri^ sin 



(5) 
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assumes the form 

„ Nabcd 
Kabcd = 



IIIabcd IIabcd 



(6) 



where the scalar quantities Nabcd, IIIabcd and IIabcd 
are of the algebraic type N, III and //, respectively, 
in the Petrov classification of the curvature tensor [llj. 
and r is the parameter distance along the principal null 
direction d/dr. Eq. © displays the peeling property [l^ 
of the curvature tensor, showing that indeed RT is the 
exterior gravitational field of a bounded source emitting 
gravitational waves. For large r we have 



Ral 



Nae 



(7) 



so that at large r the gravitational field looks like a 
gravitational wave with propagation vector d/dr. The 
nonvanishing of the scalars Nabcd is an invariant cri- 
terion for the presence of gravitational waves, and the 
asymptotic region where C'(l/r)-terms are dominant de- 
fined as the wave zone. The curvature tensor components 
in the above basis that contribute to Nabcd, namely, to 
the gravitational degrees of freedom transversal to the di- 
rection of propagation of the wave, are i?0303 = — -R0202 = 
-D{u, e, (j>)/r+0{l/r^) and i?0203 = ~B{u, 9, <?i)/r+0(l/r2) 
where 



D{u,9,<f>) 



du 



K 



1 



K. 



4,4, 



K 



and 



B{u,9,<l,) = 



K'^ sine 



du 



,KeK4 



cot 



K 



(8) 



(9) 



From (O we can see that the functions D and B contain 
all the information of the angular, and time dependence 
of the gravitational wave amplitudes in the wave zone 
once K{u, 9, (f>) is given. D and B actually correspond to 
the two independent polarization modes of the gravita- 
tional wave, transverse to its direction of propagation at 
the wave zone. 

The field equations present two stationary solutions, 
which will play a crucial role in our future discussions. 
The first is the Schwarzschild solution corresponding to 



K — Ko = const, A = I/Kq 
and mass Mschw — ttioKq. The second is 

Ko 



K{e,(b) = 



cosh 7 -I- (n ■ x) sinh 7 ' 



(10) 



(11) 



The K{9,(j)) function ([TT|) . which depends on three pa- 
rameters, is a K-transformation of t he g eneralized Bondi- 
Metzner group discussed by Sachs fisj (the BMS group) 
and represents the general form of Lorentz boosts con- 
tained in the homogeneous orthochronous Lorentz trans- 
formations of the BMS. 

The Bondi mass function of this solution is given by 
m{9,(p) = ■moK^{9,(j)). The total mass-energy of this 
gravitational configuration is given by the Bondi mass 



M 



(l/47r) 



d9 m{9, (j)) sin ^ 



mo-K'o cosh 7 — moK%l \/ 1 — 



(12) 



The interpretation of ([Tl]) as a boosted black hole is 
relative to the asymptotic Lorentz frame which is the 
rest frame of the black hole when 7 = 0. 



III. CHARACTERISTIC INITIAL DATA 

As well known the initial data problem for RT space- 
times is within the class of characteristic initial value 
formulation as opposed to the 1-1-3 formulation, accord- 
ing to the classification of Yorkfl^. For RT spacetimes 
the function K{uo, 9, (f) given in a characteristic surface 
u — corresponds to the initial data to be evolved via 
the RT equation Our task is now to exhibit an ini- 
tial K{uo, 9, (j)) that represents instantaneously the gen- 
eral collision of two Schwarzschild black holes (head-on 
collision or not), by extending a procedure outlined in 
Refs. 0, to construct initial data for RT dynamics in 
the non-axial case. 

In analogy to bispherical coordinates[16] in the 3- 
dim Cartesian plane E, let us introduce the following 
parametrization for Cartesian coordinates 

a sin9 sinh?) 



y = 



cosh rj + cos 9 sinh rj 

a sin9 sinh 77 
cosh f) -I- cos 9 sinh r) 
a 



(13) 



cosh rj + cos 9 sinh ri 



for z > and z < respectively. In the above < 77 < 
oo,O<0<7r, 0<(/)< 27r. In this parametrization, 
•q = rjQ corresponds to two spheres, one at z > and the 
other at z < 0, centered at {x = y = Q, z = ±a cosh 770) 
respectively, with radius a sinh 779. The Cartesian vector 
from a point P : [x, y, z) of E has length 



where x — (sin cos 0, sin sin 0, cos0) is the unit vector 
along an arbitrary direction x and n = (711,712,77.3) is a 



constant unit vector (satisfying n\ - 



1); also Kq 



and 7 are constants. We note that PT|) yields A = ^/Kq, 
resulting in its stationary character. This solution can 
be interpreted as a boosted black hole along the axis 
determined by the unit vector n with boost parameter 
7, or equivalently, with velocity parameter v — tanh7. 



r{ri,9) 



cosh f] — cos 9 sinh r) 
cosh rj + cos 9 sinh 77 



(14) 



For 1] = 00 the spheres degenerate into the planes z = 
and z = ±00. The usefulness of this parametrization will 
become clear in what follows. We note that the Carte- 
sian coordinates are continuous functions, with continu- 
ous derivatives, of (r],9,ip). Singularities occurring are 
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the usual singularities of a spherical coordinate system. 
For future reference let us introduce the functions 



S(^±){r],9,<p,n) — \J cosh 77 ± (n • x) sinh j] 



(15) 



where x — (sin Q cos </>, sin 9 sin 0, cos ff) and n = 
(ni, 712,^3), with + + n§ = 1. In the above 
parametrization (jl3p . the flat space line element ds^ = 



J 



(da;)^ + (dy)^ + (dz)^ is expressed as 



d77 + sinh 7; (dQ + sin^ 



(16) 



We now take S as a spacelike surface of initial data, 
with geometry defined by the line element 



ds' =a'i^'(r; + 7o,6',<^) 



d-q^ + sinh^(77 + 7o)(d6l^ + sin^ ed<3^'^) 
I 



(17) 



where 70 is an arbitrary parameter. By assuming time- 
symmetric data (namely, S a maximal slice with zero 
extrinsic curvature) we obtain that the Hamiltonian con- 
straints reduce to ('^^i? = 0. With the substitution 
iiT = the constraint equation reduces to the Laplace 
equation 



^— (-l-esine) +($'sinh'(r; 



-70, 



^0 -I- - sinh (ry 4- 70)$ 



0, 



(18) 



where a prime denotes derivative with respect to 77. It 
is not difficult to verify that the functions 



$ = 



5(±)(77 + 7o,^',<^, 



(19) 



satisfy Eq. (fTS]) and, with respect to metric (fT7|) . cor- 
respond to flat space solutions (zero curvature). It then 
follows that 



$ = 



5(„)(»7 4-7o,6',(/>,n) 



0-2 



5(+)(r? -I- 70, 6l,(^,n) 



(20) 



is a nonflat solution of (jl8l) . where a.\ and are ar- 
bitrary positive constants. The nonflat 3-dim geometry 
defined by (^0]). 



d77^ -I- sinh'^ (77 -I- 'yo){d6^ + sin^ ( 



(21) 



is asymptotically flat with a form analogous to that of 
the 3-dim spatial section of the Schwarzschild geometry 
in isotropic coordinates, as we proceed to show. 

Without loss of generality we take in (I^Hl) n — (0, 0, 1), 
that corresponds to choose the z-axis along n. In this 
instance a straightforward manipulation shows that, for 
77 >> 7o, the metric (1211) can be rewritten as 



ds^ = 02-1- 



aai ^cosh ri — cos 9 sinh 77 



''('7,^) -y/cosh 77 - (n ■ x) sinh 77 



ds 



2 

flat- 



(22) 



Now to probe the asymptotic structure of the metric 
let us consider 77 very large and, for this 77, points 
(a;, y, z) whose distance from the origin is also very large. 



namely, when (77 — ^ 00, 9 ~ tt). In this asymptotic limit, 
returning to Cartesian coordinates, the 3-geometry ([2^ 
can be given in the approximate form 



2 M(i) 



(23) 



where we fixed the scale of bispherical-type coordinates 
by taking 2^/2a — 7719(0:1 + a2)\/l + oi^Ioly. The 
Schwarzschild mass — raf^ioLx -I- a-i). 

From the above construction we can now extract ini- 
tial data for the RT dynamics, which has its initial value 
problem on null cones. Based on the initial data formu- 
lation on characteristic surfaces proposed by D'Inverno 



and Stachel[17l Il8j - in which the degrees of freedom of 



the vacuum gravitational field are contained in the con- 
formal structure of 2-spheres embedded in a 3-spacelike 
surface - we are then led to adopt the confornial struc- 
ture given by the conformal factor ([20]) defined on the 
surface 77 = 0, 



+ ■ 



02 



5(_)(7o,6i,<j!>,n) 5(+)(7o, 6*, <?!>, n) 



, (24) 



with n = (0, 0, 1), as initial data for the RT dynamics. 
This conformal structure is to be extended along null 
bicharacteristics and propagated along a timelike con- 
gruence of the spacetime via RT dynamics. A restricted 
spacetime may then be constructed locally as the prod- 
uct of the two-sphere geometry times a timelike plane 
(m, f) generated by a null vector d/df and a timelike vec- 
tor d/du with geometry dcr^ — a^{u, f, 9, ip)du'^ + 2dudr. 
The four geometry is then taken as 



ds^ = a^{u,r,9,(f>)du'^ + 2dudr 
- f'^K^(u.9.(h)(d9^ +SU 



(25) 



Eq. (l25t is the RT metric, the dynamics of which (ruled 
by Einstein's vacuum field equations) propagates the ini- 
tial data (|24l) forward in time from the characteristic ini- 
tial surface u = uq- We note that Einstein's vacuum 
equations demand that the function a^(u, f, 9, 0) has the 
form given in 

The interpretation of the asymptotically flat ini- 
tial data (|24l) as two instantaneously interacting 
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Schwarzschild black holes boosted along the z-axis is 
now discussed, based on perturbations of the RT met- 
ric ([25l) constructed with such data. As shown in Sec- 
tion II, for ai — the data correspond in (|25p to a 
static Schwarzschild black hole (with the total Bondi 
mass mo(a2)^ cosh7o) boosted along the direction de- 
fined by the unit vector n, with v = tanh7o. For ai ^ 0, 
with ai << a2, the configuration is no longer static and 
cannot therefore be a black hole, but can still be inter- 
preted as an initially perturbed boosted Schwarzschild 
black hole. Conversely the same consideration holds for 
[ai ^ 0,Q!2 = 0) and ai with a2 << ai, the latter 
case corresponding also to an initially perturbed boosted 



Schwarzschild black hole. In this sense we associate the 
perturbation with a black hole of relative small rest mass, 
boosted along the direction n (or n) in non head-on colli- 
sion with a larger black hole boosted along the direction 
fi (or n) . The initial infalling velocity of each black hole 
considered individually is given hy v — tanh 70 . 

Without loss of generality, in the remaining of the pa- 
per we fix a2 = 1 and drop subscripts 1 and of the pa- 
rameters ai and 70, respectively, in order to avoid over- 
cluttering in formulae and Figures. In this instance, the 
initial data (|24|) will in principle contain four indepen- 
dent parameters, namely, (a, 7, n) with nl + n^+n^ = 1, 
and assumes the form 



J 



1 



+ ■ 



V cosh 7 + cos 61 sinh 7 ^cosh 7 — (n • x) sinh 7 
I 



(26) 



We must comment that, since only two black holes are 
involved in the collision, the data ((26|) should actually 
depend (besides a and 7) only on the incidence angle 
Pq between the two initial black holes, determined by the 
scalar product nz -n. Therefore we should expect that (i) 
the resulting dynamics of the system would not be altered 
under a rigid rotation of the two initial black holes about 
the z-axis, which could be used to locate n of the initial 
data in the x — z plane; and consequently (ii) the initial 
data will result in a planar dynamics of the collision, a 
fact that would save a lot of computational effort. The 
above statements (i)-(ii) will indeed be confirmed by the 
numerical results obtained from the evolution of the data 
(PS)) . as shown in section VI, and establish the planar 
nature of a general non head-on collision of two black 
holes. 

Finally we should remark that, in the full Bondi- Sachs 
problem, the analysis of field equations in the 2 + 2 
formulation^, Il7| shows that the two news functions 

ci^\u, 9, (j)) and c^u\u, 0, 0) are part of the initial data to 
be prescribed for the evolution of the system 19] . How- 
ever for the RT dynamics the news are already speci- 
fied once the initial data for the RT equation, namely 
K{uq, 9, (j)), is given and consequently K(u, 9, (p) is given 
for all u > uq from the numerical evolution of the data. 

The initial data (j26p will be evolved numerically - via 
RT dynamics - up to a final configuration that corre- 
sponds to a remnant Schwarzschild black hole boosted 
along the axis determined by the unit vector nj and 
having the form of the solution (|lip . as we discuss in 
the next Section. A numerical code using the Galerkin 
spectral method is was implemented to integrate the non- 
axisymmetric RT equation (U), the basis of which is de- 
scribed in the next Section. 



IV. NUMERICAL EVOLUTION OF THE DATA 

We proceed now to discuss the numerical evolution of 
the initial data (j26l) via the non-axisymmetric RT equa- 
tion ([4]). Throughout the present Section the variable 
9 will be expressed in terms of the variable x = cos 9. 
The numerical integration of the non-axisymmetric RT 
equation is performed using a Galerkin spectral method 
which is now described in detail. In the present pro- 
cedure we rewrite the equations ([3])- ([4]) using the vari- 
able P{u,x,(j)) = 1/K{u,x,(j)) (instead of the former 
K{u,x,(f>)), an approach already adopted in Ref. [g^. 
We obtain that 



A(u, X, ct>) = {l- x^) [PP^, - P|] - 2xPP,. +P^ + 



+ 



^ 'PP^4, - Pi] 



(1-2:2) 



and 



P(u, X, 



12mo 



(1 — X )\xx — 2x\x + 



^4,4, 



(1 - a;2 



(27) 



(28) 



where a dot and the the subscripts x and denote, 
respectively, derivatives with respect to u, x and 4>. 

The Galerkin method establishes that P(u, x, 0) can 
be expanded in a convenient set of basis functions of a 
projection space by which we can reduce the RT partial 
differential equation ([25]) into a finite set of nonlinear cou- 
pled ordinary differential equations. This set of equations 
constitutes an autonomous dynamical system, the dimen- 
sion of which depends directly on the truncation in the 
Galerkin method to approximate the RT dynamics. As 
the conformal function P{u, x, (f) defined on the 2-sphere 
is assumed to be sufficiently smooth, we can use the real 
Spherical Harmonics (SHs) [llj as the appropriate ba- 
sis that better approximates our desired solution. Before 
starting with the numerical scheme, we first present some 
definitions and properties of these functions. 
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Unlike the complex SHs, the real SHs consist of a break 
of the complex SHs into their sine and cosine parts as 
follows, 

(s,0iZ,m)(+) = ^+^Yr{x,c^) = Wi,mPr{x)cos{m<j,) 

= <+){/, m|a;,<^), 
{x,,^|;,m)(-> = <-'yr(a;,0) = W^i,mPr(a;)sin(m<j!.) 

= '■-^ {l,m\x,^), (29) 

where Pl^{x) are the associated Legendre functions de- 
fined by 



p 



and the normalization factors Wi,m are given by 



1 {21 + I) {I - my. 

2 HTTmy. 



(30) 



(31) 



From these formulas we see that the case m = gives 
us the basis functions for the axial case, the Legendre 
polynomials Pi{x) = {x\l) plj . 

The real SH orthogonality relations are given by 



J -1 
1lT /■! 



(«, mjx, (^) (a;, 01/', m') '±>da; 



(±¥r(a;,0)(±¥,y^'(a;,0)dx 



/ J - 1 



(32) 



and 



^27r /.I 



J J ^^'> {l,m\x,4>) {x,<l}\l',m')^^^dx . 



{±) 



Yr{x, (t>)''^^YiV {x, 4>)dx d<i> = 0, (33) 



and are fundamental in the treatment of the Galerkin 
method. Here the following completeness relation is also 
used, 



I I \x,4>){x,4>\ dx d(l) = 1. 



(34) 



Now we are able to construct our numerical Galerkin 
scheme, that will allow us to obtain an accurate approxi- 
mated numerical solution for (j28p corresponding to given 
initial data. As every twice continuously differentiable, 
suitably periodic real function defined on the surface of 
a sphere admits an absolutely convergent expansion in 
terms of the SHs, let us consider the expansion 

JVp 



P(„)(«, x,^) = Y. \aUu)^+^Yi\x, cj>) + 

JVp n I 



I 



+ J2 Bi,^{u)<--^Yr{a 



(35) 



as our approximated solution for P{u, x, (f). Here Np is 
a positive integer that defines the truncation order of the 
Galerkin method. By using the orthogonality relations 
([5^ and ([55)1 . the modal coefficients Ai^miu) and Bi miu) 
are given by 

ALmiu) - J J ''*^Yi"'{x,(t>)P{a){u,X,(l>)dxd(l> 

= (+)(;,m|P(„)), (36) 



Bi^,„{u) = 1 / ^~^Yl^{x,(j>)P(a){u,x,(j>)dxa 
- (-) 



(/,m|P(„,) , 



(37) 



Inserting the expanded solution ([35]) in the RT partial 
differential equation (1281) we obtain an autonomous dy- 
namical system of dimension {Np + 1)^ for the modal 
coefficients {Ai,jn{u), Bim{u)), namely, 

Al.rniu) = Al^jn[Ai^jn{u),Bi,jn{u)^, 
Bl,ra{u) = Bl^,n(^Ai^rniu), Bi„iiu)j , (38) 

where Ai^m and Bi^m are polynomials of order {Np+1)'^ in 
Ai^miu) and Bi „i{u) and a dot denotes here u-derivative. 
The Galerkin scheme guarantees that the projections 
of {P{u, X, (j)) — P(a) (?i, .X, (/))) onto each basis function, 
namely, < {P{u,x,<l)) - P(a){u,x,(j))),Y/^{x,(t)) > ap- 
proach zero when Np — > oo so that ([35)) approaches 
the exact solution P{u, x, (j)), in the sense of the norm 
of the projection basis space of the SHs. The initial 
conditions {Ai^miuo), Bim{uo)) to be used to integrate 
the dynamical system (|38[) are provided by the initial 
data K{uo,x,(p) constructed in last Section (Eqs. (1^ ). 
These initial values are obtained from the Galerkin de- 
composition of P{uo,x,(j)) (cf. Eq. ([35)) ) evaluated ac- 
cording to ([5S)) . 

As a matter of fact, the projections that transform the 
RT equation into the dynamical system ()38p demands ac- 
tually that we replace P{u,x,4)) by P(a)iu,x,(j)) in ([77)) 
and ([^5)) . But in doing so we excessively increase the 
computational demand, overloading the computer data 
storage. To circumvent this problem, we divide the RT 
equation into two parts that we called constraints: the 
A-constraint and the P'^-constraint. The process will con- 
sist then in also expanding - in the same way as in ()35p 
- the expressions for A and P^, 

A(„)(u, x,<P)^Yl ^'^M^''^Y°{x, 4>) + 

1=0 
Nx I 

i = l m = l 
Nx I 



iVp3 

P\^^{u,x,(f,) = l^Ci,o{u)'^+^Y['{x,^) + 
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1 = 1 7n = l 



Here, TVa and Nps are integers that define the trunca- 
tion orders for the respective expansions. Throughout 
the paper, we take Np = N\ = Nps, but these orders 
are totally independent and can be taken with diferent 
values. Now, we have a new two set of modal coefficients 
{{Cim, A,m} and {Ei^jri, Fi^rn}) that wiU be uniquely de- 
termined by the first employed set of modal coefficients 
{Ai^m,Bi,„i} by 



* + ' {lMP\a)) 



Eum{{Al^m{u), Bl^m{u)}), 

(;,m|A(„)) 

FLm{{ALm{u),Bl,^{u)}). 



(41) 



The RT dynamical equation will have then a new form: 
the left side will be written by the it-derivative of P(a) 

Np 

P(,)(«, X, Cl>)teft = ^ -A,o{u)^ + ^Yt°{x, 0) + 
1=0 
Np I 



1 = 1 m = l 
Np I 



(42) 



and, its right side, will be given by the replacement of 
A and in ([25]) by the constraint expressions ([M]) and 

601). 



P(„)(U, X, 4>)right 



(ay) 
V2mo 



X 



Mi 



Finally we project both sides to get the relation be- 
tween the M-variation of the set {^(^^(m), 5i,m(w)} and 
the new modal coefficients 

Ai,^{u) = m!P(„)) 

= Al,m{{Cl,m{u),Dl^rn{u),El^rn{u),Fl^rn{u)}) 

Bi,r^{u) = <^Z,m|P(„)^ 

= Bl^m{{Cl^m{u), P;,m(u), Pi,m(^t), Pi,m('i)}l43) 

This completes our scheme to obtain the dynamical sys- 
tem version of the non-axisymmetric RT equation (|28p. 
To solve it, we use a fourth-order Runge-Kutta recur- 
sive method adapted to our constraints. So we have 
to get the initial values for the modal coefhcients, say 



{Ai,™(Mo),Si,m(uo)} 



Ai,m{uo) = (;,m|Po(„)) , 



(44) 



The effect of the truncation Np on the initial data 
may be evaluated by the relative error TZ{uq) between 
the approximate and exact expressions given by pSI) and 
P{uo,x,phi) = l/K{uo,x,phi) (cf. (|26|)). 



TZ{uo,x) 



P{uo,X,(t>) - Pi^a){uQ,X,( 
P{uo,X,(j)) 



(45) 



In all the numerical experiments of the paper we have 
adopted Np = 7; with this truncation the relative error 
above is of the order of, or smaller than 10~*, for all —1 < 
X < 1 and < (/) < 27r. The RT dynamics furnishes us 
with a further test to check the accuracy and reliability of 
our numerical codes. In fact, for any sufficiently smooth 
K{u, x) we have that the quantity 



C(«) 



P {u, X, (j)) dx, 



(46) 



is conserved along the dynamics, namely, duCiu) — 0. 
Evaluating its exact value from the initial data (|26p and 
at distinct steps of the computation we that ICCwo) — 
C{u)\ < 10"^*^ for all our numerical experiments and for 
all the sampled values of m > wq, for the adopted trunca- 
tion Np = 7. 

Our numerical experiments are realized with the initial 
data (|26p - which corresponds to one initial black hole 
boosted along the 2-axis and the second black hole mov- 
ing along the direction determined by n = (711,712,713) 
with n\+n\+n\ — 1- having in principle four indepen- 
dent parameters. 

We vary a in the interval (0, 1.5) for 7 = 0.2 and 
several values of ni. Exhaustive numerical experiments 
show that for a sufficiently large computation time Uf 
all modal coefficients become constant, namely, at Uf 
the emission of gravitational waves is considered to have 
ceased. For Uf we actually have that \Ai^rniuf + h) — 

Amiuf)] < 10-1", \Bi^miuf + h)- Bi^^(uf)\ < 10-12 

for all I = 0...7, where h is the integration step. 

From these modal coefficients Ai^rn{uf) and 
Bi,rn{uf) we reconstruct the final configuration 
P{uf,x,(/)):^P(^a)iuf,x,(l>) that, in all cases, can 
be expressed as 



P{Uf,X,(j)) = 



^o,o(w/) Aifijuf) ^^^^ 



2 ' 2 
+ y4i,i(u/) sin S cos 

+ Pi, i(u/) sin 6* cos + (10-1°) 

- (^"^'i 7/ + n/ • X sinh 7/ ) . (47) 



The rms error of the second equality in (|47)) is of the 
order of, or smaller that lO'i^. The expression in the 
second equality of (|T7)) corresponds to a boosted black 
hole along the direction = (niy, 712/, 713/) with boost 
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parameter 7/ and rest mass irioKj, a confirmation of the 
Chrusciel- Singleton theorem p^]. 
Prom ((47|) we can identify 



2 cosh 7/ 
2 723/ sinh7/ 
nif sinh7/ 
n2f sinh7/ 



(48) 



K 



f 



up to 0(10^^^), from which we can read the final param- 
eters of the remnant stationary black hole {Kf,jf,nf). 
It results 



nif = 



2Ai,iK) 



(^1 + (Ml:iM)2 + (^^liM) 

Ai,o{uf) Aifiiuf) 



2Bi,i(^t/) 



n3f 



-1/2 



7/ = tanh (— - — ) 
n-if Aofi 

2 

Kf = - — cosh 7/. 

^0,0 



(49) 



(50) 



An alternative to evaluate the parameter Kf is the use 
of the initial data in the conserved quantity (|46p , namely, 



1 /■2ir /■! 
i / , , / ,,2/ 



if/ = ( — I d<j) I K {uo,x,(j)) dx 



1/2 



(51) 



The agreement is within 10 significant decimal digits. 
This is also a test of the accuracy of (ni/)^ + (n2/)^ + 

One of the basic results to be extracted from our nu- 
merical experiments are the values of (Kf,^f,n^f), for 
each of the independent parameters (a, 7, n) of the initial 
data. The former are the basic parameters of the remnant 
that - together with the initial data of the system and 
the function K[u^d^(j)) for all uq < u < Uf~ allow us to 
evaluate quantities which are characteristic of the radia- 
tive transfer processes involved in the gravitational wave 
emission. Finally we remark that the integration of the 
dynamical system (H51) ) used a fourth-order Runge-Kutta 
recursive method (adapted to our constraints) together 
with a C-I-+ integrator. Unless otherwise stated, all our 
numerical results are restricted to the choice 7 = 0.2, 
corresponding to the initial infalling collision velocity 
V ~ 0.1973. 



V. THE PLANAR NATURE OF THE GENERAL 
NON-HEAD-ON COLIISION 

From the numerical experiments done to check the long 
time evolution of the initial data and to obtain the 
basic parameters of the remnant black hole, we observe 
that the final direction of remnant velocity is contained 
in the plane determined by the the two directions of the 
initial individual colliding black holes. For illustration 
let us consider the evolution of the initial data (|26p , with 
Uz = (0, 0, 1) corresponding to the direction of the initial 
colliding black hole (boosted along the positive 2;-axis), 
with n = (0.26,0.07,0.963068) corresponding to the di- 
rection of the second initial colliding black hole and the 
several values of a = 0.1, 0.1, 0.20, 0.3, 0.50. The result- 
ing final direction of the remnant is given by (cf. (|49|) ) 



'^ 3.670544 10~^ 9.999003 10"^) 
"^ 8.251592 10"^ 9.994961 10"^) 
,0.000000 10"^ 9.000000 10"^) 



n/(Q = 0.1) ~ (1.363345 10" 

n/(a = 0.2) ~ (3.064877 10" 

n/(a = 0.3) ~ (0.000000 10"^ 

n/(Q = 0.5) ~ (1.197374 10"\ 3.223699 10"^ 9.922820 10"^). 



In all cases we obtain 



• n A n/ ~ 0. 



(52) 



Within machine precision, (|52p is valid up to 0(10^^"), 
and it also holds for all numerical experiments done in 
the paper considering a large domain of the initial data 
parameters, implying that the dynamics is indeed planar, 
as discussed at the end of Section III. The result ([5^ 
could also be considered as an additional test for the 
accuracy of our numerical code. 

This fundamental result - which will be used to save 
computational efforts in the dynamical system evolution 
of the RT non-axisymmetric dynamics - has an impor- 
tant physical counterpart connected to the fact that the 
linear momentum flux carried out by gravitational waves 
is confined to the plane of the initial collision. In other 
words, the linear momentum of the merged system in the 
direction nzAn is conserved, as discussed in Section VIII. 

Let us then consider a rigid rotation of the two initial 
black hole system about the direction of motion of the 
first initial black, namely, the z-axis. This rotation will 
be chosen so that the plane of the initial collision, deter- 
mined by the two directions = (0, 0, 1) and n, coincide 
with the X — z plane. In this case the direction of motion 
of the second initial colliding black hole is now expressed 
by the unit vector n = (rii, 0, ^2), with n\ -\- n\ = 1. In 
this instance the initial data (|26|) can assume the simpler 
form 



V cosh 7 -I- cos 6 sinh 7 ^ cosh 7 — (cos po cos 6 + sin po sin 6 cos (j)) sinh 7 

I 



(53) 



This initial data contains now three independent pa- rameters (a, 7, po) only, where po is the angle formed by n 
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and the z-axis and is used to parametrize ni = cos po and 
= sin pq . Again, as already established, the motion of 
the resulting remnant black hole will then be restricted 
to the X — z plane for any values of the initial data pa- 
rameters. From the preceding discussions of the present 
Section we have that the initial data ((53|) describes a gen- 
eral non head-on collision of two black holes in the realm 
of RT dynamics. 

Now, for exhaustive numerical tests with the initial 
data ([55]) and a large range of parameters (a,po) we 
verify that: (i) the initial modal coefficients Bi^m{uo) < 
10^^^, and (ii) its evolution via the RT dynamical system 
(|43| maintains -Bz,m(u) ^ 10~^^, for all {l,rn) in their al- 
lowed range, and for all mq < u < Uf. This corresponds - 
within the precision of our computation - to Bi^m{u) = 
for all u. Obviously, from the second equation (P^ . we 
have that the unit vector determining the direction of 
the remnant will always have the form n = (ni/,0,n3/). 
Therefore our computational task can be simplified by 
taking all the Bi „i coefficients equal to zero, which is 
equivalent to restrict all our expansions in the Galerkin 
method decomposition to the cosine series only. In the 
remaining of the paper we have adopted this procedure 
in the numerical evaluations of RT dynamics for the data 
([55)1 . We have however made sample tests by verifying 
the absolute differences in the physical results originat- 
ing from data either from the complete decomposition or 
the decomposition using the cosine series only. The dif- 
ference remains always of the order of, or smaller than 
10^1^. The angle po (0 < po < 7r/2) will be denoted the 
collision angle, the limiting cases po — corresponding 
to a head-on collision and po = 7r/2 corresponding to a 
right angle collision. 

It will be possible to follow the full evolution of ([55)1 
by considering from very small a up to those values for 
which the nonlinearities start to play an important role 
in the dynamics. We will restrict ourselves to the range 
< a < 1.2 The range of the initial incidence angle in the 
numerical experiments will be < po ^ 90°, the limits 
corresponding respectively to a head-on collision and an 
orthogonal collision. 



VI. THE BONDI-SACHS FOUR MOMENTUM 
AND CONSERVATION LAWS 



satisfy the Bondi-Sachs boundary conditions. It should 
be noticed that in the RT coordinate system the pres- 
ence of the term 2rKu/K does not fulfill the appropriate 
boundary conditions. Although the coordinate transfor- 
mations from RT coordinates to Bondi-Sachs coordinates 
cannot be expressible in a closed form (they are given by 
an infinite series in powers of r~^) p^ . their asymptotic 
expansion allows us to obtain the form of the required 
physical quantities. In this section we will restrict our- 
selves to the Bondi-Sachs energy-momentum for the RT 
spacetimes as well as its conservation laws. Our deriva- 
tion for the non-axisymmetric case[l^ follows closely the 
work of Gona and Kramer [2^ done for the axisymmetric 
case. 

From the supplementary vacuum Einstein equations 
Ruu = 0, RuQ^o, and Ru<s>=o in the 2-1-2 Bondi-Sachs 
formulation [3, lg| (where (U, R, 9, $) are the Bondi-Sachs 
coordinates), we obtain in RT coordinates 



dm{u, 6, (j)) 



+ 4c 



du 

(2) COS( 



= -A'(c« +c(^' 



2\ ld_ 

^ 2 9^ 



3ci^' cot ( 



ZC + Cgg -I- Cg . 



(54) 



where to(u,0, is the Bondi mass function and 

c'ii\u,0,(j)) and ou\u,0,(l)) are the two news functions 
for the non-axisymmetric case[l9|. corresponding to the 
two modes of polarization of the gravitational waves. The 
extra factor K in the first term of the second-hand-side 
of Eq. ([5l)) comes from the transformation 



lirrir 



' du 



1 



(55) 



U being the Bondi time coordinate. For the news satisfy- 
ing the appropriate boundary conditions c^^^ = c*^^-* = 
and Cg^-* = Cg^"* 

Sachs four-momentum conservation 



Cg ' = at = 0,7r , we obtain the Bondi- 



du 



- 



(56) 



where the Bondi-Sachs four-momentum P^(u) is defined 



as 



and 



1 p277 p77 

^''(u) = — d4> m(u,e,(j)) F sine de, (57) 

47r Jo Jo 



RT spacetimes describe the asymptotically flat exterior 
gravitational field of a bounded system radiating gravita- 
tional waves and in this sense they are in the realm of the 
2-1-2 Bondi-Sachs formulation of gravitational waves in 
General Relativity [Vl IS]- Furthermore initial data for RT 
dynamics are prescribed on null characteristic surfaces. 
Therefore suitable expressions for the physical quantities 
to be used in the description of gravitational wave emis- 
sion processes and its conservation laws must be derived. 
To exhibit such expressions it is necessary to perform a 
coordinate transformation from RT coordinates used in 
([T]) to a coordinate system in which the metric coefficients 



d(l> I K ci^>^ + c$f'^ ) sine de, (58) 



corresponds to the net fiux of energy-momentum carried 
out by the gravitational waves emitted. In the above 
l^^ = (1, — sin6'coS(/), — sin^sini/), — cos 6*) is a null vector 
relative to an asymptotic Lorentz frame at infinity. We 
note that the last term in ()54|) vanishes in the integrations 
due to the boundary conditions of the news. 

For p — Eq. ([56]) yields the Bondi mass formula 



dMsiu) 



-Pw{u) 



(59) 
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where Mb{u) is the Bondi mass at a time u and 



d(f) 



K 



,(1) 



(60) 



is the power extracted from the system by the gravita- 
tional wave emission in a time u. The total energy E\y 
carried out of the system by the gravitational waves is 
given by the time integral of (j60p up to m/, and cor- 
responds to the total Bondi mass extracted from the 
system, Mb(wo) — Mb(m/) = Ew^ where Msiuf) — 
nioKj cosh 7/ is the Bondi mass of the remnant. 



For fi — x,y, z, Eqs. (|56p determine the liner momen- 
tum conservation of the system, 



d F{u) 

du 



Pw(w) 



where the vector 



'w[u) 



47r 



•-Ti (--7/ 



(61) 



sinOde. (62) 



out by 
above x 



the 



is the net momentum flux carried 
gravitational waves emitted. In the 
(sin cos 0, sin 9 sin cos 9). 

We are now able to analyze the radiative processes that 
lead the merged system from its initial configuration to 
the final configuration of the remnant black hole. 




0.00004 



0.00003 



0.00002 



0.00001 



u/mo 

FIG. 1: Linear-log plot of the energy flux (power) carried out 
by gravitational waves as a function of u/mo, for a = 0.1 
and 7 = 0.2 and several values of the collision angle po. We 
see that, for this relatively small value of a, the gw emis- 
sion corresponds to a pulse of short duration. Aw/mo ~ 3.5, 
and its initial intensity decreases as po increases. The case 
Po = (head-on collision) constitutes an upper bound, while 
the orthogonal collision (po = 90°) constitutes a lower bound 
for the total energy emitted (the area below the curve), in 
accordance with the efficiency behavior. 



VII. THE POWER EXTRACTED BY THE 
GRAVITATIONAL WAVE EMISSION AND THE 

EFFICIENCY OF THE MASS-ENERGY 
RADIATIVE TRANSFER: UPPER AND LOWER 
BOUNDS 

We now examine the energy extracted from the merged 
system by the gravitational waves. From our discussion 
in Section VI on the Bondi-Sachs conservation laws, and 
specifically from Eqs. (|59|) and (|60|) we have that the 
power emitted by the system in a time u is given by 

Pw(u) = dEwiu) / du 

= ^ dcf) 1^ Jf^ci^^Vc'f^^^ sinSde. (63) 

In Figs. [Hand [2] we plot the power emitted Pw{u) as a 
function of u/rriQ for a = 0.1 and a — 0.6, respectively, 
and 7 = 0.2 fixed, and for several values of the incidence 
angle po- The total energy emission Ewiuf) in each case 
is measured by the area below the respective curve ~ 
in accordance with (j60l) . Two important features are 
to be noted in the Figures. First, we can see that a 
head-on collision (po = 0) constitutes an upper-bound 
for the total energy emitted, as well as the orthogonal 
collision (po = 90°) constitutes a lower-bound for this 
energy. This pattern is typical for any < a < 1 and for 
all 7. Actually in our numerical experiments we verified 
this important result for < a < 1, consistent with 
the curves of the efficiency of mass-energy extraction by 



gravitational waves, as we will discuss later in the present 
Section. 

Second, although the curves present the same behavior 
for distinct a's, we can see that for a small (cf. Fig. [T] 
for a = 0.1) the emission corresponds to a short pulse of 
gravitational waves. We indeed obtain that, for any of 
the incident angles po (here including the head-on colli- 
sion) , the initial power emitted decreases by three orders 
of magnitude in an interval of time Au/mg ~ 3.5. On 
the contrary, for a = 0.6 (cf. Fig. ^ the initial power 
emitted (which is one order of magnitude smaller than 
in the case a — 0.1) decreases by three orders of mag- 
nitude in an much larger interval Au/mg 72 for all 
the incidence angles po considered. In fact, analogous to 
the case of head-on collision analyzed in Ref . [2] , there 
is a threshold value of a ^ 0.67 that separates regimes 
of short bursts of gravitational waves from a regime of 
quiescent long time emission. 

For a in the interval (0, 1.5] and for several values of po 
in the interval (0,7r/2) we have determined {Kf,"ff,pf) 
which characterize the final boosted black hole configu- 
ration. From these parameters of the black hole remnant 
and the total radiated energy Ew(uf), the efficiency A 
of mass-energy extraction by gravitational wave emission 
can be evaluated. According to Eardlev[2^ the efficiency 
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0.00015 



^0.00010 
ft, 



FIG. 2: Linear-log plot of the energy flux (power) carried out 
by gravitational waves as a function of it/mo, for a — 0.6 and 
7 = 0.2 and several vaues of the collision angle pQ. Analogous 
to the case of a = 0.1, the head-on collision constitutes an 
upper bound, and the orthogonal collision a lower bound for 
the total energy emitted in accordance with the efficiency be- 
havior. The total energy emitted is spread on a larger interval 
of time, Au/mo ~ 72. 
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1e-06 
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FIG. 3: Log-linear plot of the efficiency A for the range 
ai = (0, 1] and several values of the incidence angle po. The 
initial parameter 7 — 0.2 is fixed but the distribution of 
points for other values of 7 exhibit an analogous behavior. 
The points are connected for a better visualization. Values 
of A are redundant due to the relation A (a) = A(l/a) for 
a € (0,00). 



A is defined as 



Mb(0) - Meiuf) 
Mfl(O) 



(64) 



where Mb{uo) is the initial Bondi mass and Msiuf) is 
the final Bondi mass. We can alternatively express A as 



A = 



Ew{uf) 



moK'j cosh7/ -I- Ew{uf) 



(65) 



In Fig.l we display the log-linear plot of the efficiency 
A versus the mass ratio parameter a for 7 = 0.2 and 
several values of the collision angle po varying between 
Po = (head-on collision) and po — 90° (orthogonal col- 
lision). We see that the case of head-on collisions consti- 
tute an upper bound for the efficiency A while the case of 
orthogonal collisions constitute a lower bound for A, in 
accordance with the power extraction behavior discussed 
above and illustrated in Figs. [3] and SI The efficiencies 
also decrease as the initial angle of collison po increases. 

We can also derive that, due to the form of the initial 
data (j53p . the efficiency A has a maximum for a = 1. 
In fact, let us consider a > 1 in ([55]) and factorize a 
so that the initial data will turn out to be the prod- 
uct of times a new K(uo, 0, </>), the latter correspond- 
ing to the initial data of two colliding black holes, one 
with mass parameter (1/a) boosted along the negative 
z axis and the other with mass parameters 1 boosted 
along the direction n = (sin poj 0, cospo)- By a rigid 



counter clockwise rotation of the x ^ z plane, about 
the Cartesian axis y, of an angle n — po the unity vec- 
tors n and n = (0, 0, 1) will transform respectively into 
n ~ (0, 0, —1) and n = (— sinpo, 0, — cos po) reproducing 
([53)1 for a — )■ l/a. Now it is not difficult to see that, un- 
der the rescale of K{uo, 6, (p) by a factor a^, both Msiuf) 
and Ew{uf) rescale with so that the efficiency A sat- 
isfies 



A(a) = A(l/a) 



(66) 



for a G (0,00). The rescale of Ew with a® uses Eq. 
and also that du rescales with , the latter derived 
directly from the RT equation. 

Since we can see numerically that A increases with 
a in the range a = (0, 1) it necessarily follows that A 
has a maximum at a = 1 for any < po ^ 90°, as 
illustrated in Fig. 01 This property holds also for any 
7. Also we have that the difference of the efficiencies 
between the limit configurations, po — 0° (upper bound) 
and Po — 90° (lower bound), increase as a increases, 
reaching a maximum at a = 1. For the present case 
7 — 0.2 we obtain numerically that Amax(/Oo = 0°) — 
1.191990 X 10-4, Amax(po = 21°) ~ 1.115352 x 10"", 
A„,ax(po = 45°) ~ 0.872631 x lO^^, A„,ax(po = 90°) ~ 
3.029724 X 10"^. 
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FIG. 4: Plot of the efficiency A in the range ai — (0, 1.55) ex- 
hibiting the absolute maxima at a = 1 for po = 90°, 45°, 21° 
and 0°. The points are connected for a better visualizaton. 
This behavior of maxima at a = 1 actually holds for any 
< po < 90° and any 7. 



VIII. GRAVITATIONAL WAVE RECOIL: THE 
MOMENTUM FLUX CARRIED OUT BY 
GRAVITATIONAL WAVES AND THE 
MOMENTUM OF THE REMNANT 

The analysis of the linear momentum extraction from 
the merged system by the gravitational waves is made 
using the Bondi-Sachs conservation laws (1561) for fi — 
x,y, z, which reads 

Pm'(") = ^^ '^'^^ A'x^ci^'Vc'f'^^ sinSde, (67) 

where Fwiu) is the momentum flux carried out by the 
gravitational waves emitted. In accordance with the re- 
sults and discussions of Section VI the linear momentum 
of the merged system orthogonal to the plane of collision 
is conserved, which in our choice of data (I53p implies that 
~ namely, the net flux of linear momentum car- 
ried out by gravitational waves is restricted to the plane 
of collision. Therefore from (|57)) we have 

Pw{u) = d.'i^ sin^6lcos(;!> Js'^ci^'Vc'f' ^^de, (68) 

Pw{u) = '^^ cos e sin e k(^c[1'''^ + c[f'''^^d9, (69) 

as the relevant momentum fluxes contributing in the 
conservation law (|56p and determining the total impulse 
imparted on the merged system by gravitational wave 
emission. 

The total momentum conservation and the issue of 
kicks on the merged system were already analyzed in Q 
for the case of a head-on collision (po =0°). For values 
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FIG. 5: Linear-log plots of the net momentum fluxes 
and carried out by the gravitational waves, for a — 0.1, 
7 = 0.2, and incidence angles po = 15° (top) and po — 90° 
(bottom). Due to the planar nature of the collision, the mo- 
mentum flux P^ is conserved. For the case of the relatively 
small incidence angle po ~ 15° (top) we have an initial regime 
with positive P^(u), < u < Wfc, and a consequent increase 
of the linear momentum z-component of the merged system, 
dP^ {u) / du > (cf. (|67p ). Along the x axis the merged system 
is decelerated for all u until the flnal configuration is reached. 
For a large value of po = 90° (bottom) the pattern changes 
and the system is decelerated along both directions x and z. 
The initial phase with P^{u) is not present for po > 61°, as 
shown in Fig. [B] 



of the incidence angle po sufRciently small the picture is 
analogous to that of a head-on collision, with P^{u) ~ 
for all mq < u < Uf. However as po increases a negative 
net momentum flux appears along the a;-direction lead- 
ing to a decrease of the momentum of the merged system 
in that direction for all < m < u/. This is illustrated 
in Fig. [5] (top) where the net momentum fluxes P^/ and 
carried out by the gravitational waves are plotted as 
a function of u, for mass ratio parameter a = 0.1 and 
incidence angle po — 15°. We see that we have an initial 
phase where the momentum of the merged system along 
the z axis increases, due to a positive net momentum 
flux P^r{u) for uq < u < Uk, where Uk/rno ~ 0.085 is 
the time when P^ changes sign. Therefore we have a 
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dominant deceleration regime of the merged system un- 
til the remnant black hole configuration is reached. For 
a = 0.1 fixed, this basic pattern is maintained up to 
Pq ~ 61°, beyond which both net gravitational wave mo- 
mentum fiuxes, P^/{u) and P^{u), are negative for the 
whole domain of evolution uq < u < Uf. Actually for 
Pq ^ 61° the initial positive net momentum flux along z 
disappears, with the system being decelerated along both 
directions x and z in the whole domain of the dynamical 
evolution. This is illustrated in Fig. [5] (bottom) for the 
limiting case po — 90° (orthogonal collision). The behav- 
ior of the threshold conflguration po ~ 61° is illustrated 
in Fig. |6l In a forthcoming paper we show that a basic 
role of the initial phase of positive net momentum flux 
(u) is to induce a large inspiral branch in the motion 
of the center-of-mass of the merged system. 
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FIG. 6: Linear-log plots of the net momentum fluxes and 
P^ carried out by the gravitational waves, for a — 0.1, and 
the initial incidence angle po = 61° which corresponds to a 
threshold configuration where the initial phase of positive net 
momentum flux along the z axis not present, with P^{uo) = 
0. In this case we can see that the total impulse imparted to 
the merged system along the x direction is considerably larger 
than the impulse in the z direction. 



The final momentum of the remnant can be obtained 
by evaluating (|57p at u = Uf. Taking into account that 
the Bondi mass aspect for the final black hole remnant 
is given by 

mo K^t 

m u, e, 0) = — — — . ■ ,3 , 70) 

(cosn7/ + (n/ • xjsmh7/)'^ 

(cf. Eqs. and (H7)) ). we obtain by a straightforward 
integration that 

Pf^n.fPf, Pf=n,fPf, (71) 

with modulus Pf = mg Kj sinh7/. As expected we have 
py 

f . r-, 

In Fig. ujwe plot the distribution of the final momen- 
tum modulus Pf of the remnant black hole as a function 
of a, for several values of the incidence angle ranging 



FIG. 7: Plot of the modulus Pf of the momentum of the 
remnant black hole as a function of the mass ratio parameter 
a, for several values of the incidence angle of collision po- 
Points are connected for a better visualization. Values of 
Pf for Q > 1 are redundant due to the relation Pf{a) = 
a'^Pfil/a). 



from Po = (head-on collision) to po — 90° (orthogo- 
nal colhsion). We see that the head-on case {po = 0°) 
is a lower bound for the remnant moment distribution. 
This is the only case in which the modulus of the final 
momentum of the remnant becomes zero, for a = 1, as 
discussed in Ref. [Sj] where the axisymmetric dynamics 
was analyzed, making explicit the constrained character 
of a head-on collision. For sufficiently small values of the 
incidence angle pq we see that the behavior follows that 
of the head-on case, with the zero momentum at a = 1 
being replaced by a minimum at a < 1. This minimum 
turns into a inflexion point for pQ ~ 15° and disappears 
for larger values of po- Due to the dominant decelera- 
tion regime of the merged system, consequent of the lin- 
ear momentum flux Pw extracted by the gravitational 
waves, the total momentum of the system decreases to- 
wards the distribution of Pf given in Fig. [71 Finally, 
through an analogous analysis to the case of the efficiency 
A, we have that the final momentum Pf satisfies the re- 
lation Pf{a) = a^Pf{l/a) for a G (0, oo). Therefore in 
Fig. [7] the values of Pf for a > 1 are redundant. 

A comment is in order now. An important new feature 
appears in non head-on collisions connected to the final 
momentum for equal mass black holes {a = 1). Contrary 
to the case of a head-on collision and of the merging of 
black hole binary inspirals, for a non head-on collision 
of two equal mass black holes the net gravitational wave 
flux emitted is nonzero. This is the reason why the final 
momenta in Fig. [7]for a = 1 is nonzero, even if measured 
in an inertial frame with velocity Vi„ = P{uq)/Mb{uo) 
(namely, a zero-initial-Bondi-momentum frame) relative 
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FIG. 8: Linear- log plot of the net momentum fluxes P^{u) 
and P^{u) (top), and the associated total impulses Iwi''^) 
and /^(ii) (bottom), for the case of equal-mass initial collid- 
ing black holes (a = 1) and incidence angle po = 15°. The 
associated kick for this case is Vk ~ 0.9km/s, this small value 
being due to the value 7 — 0.2 adopted in our computation 
(cf. text). 



to a rest inertial frame at infinity. This is illustrated 
in Figs. [5] , where we display the net momentum flux 
Pw (u) and the associated impulse Iw (u) for the case of 
equal-mass initial colliding black holes in a non head-on 
collision with incidence angle po = 15°. 

The analysis of the kick velocities generated in these 
processes of momentum extraction (and the associated 
gravitational wave recoil) will be the subject of a future 
paper. In a related analysis (cf. Ref. 3) for head-on 
collisions we have shown that the distribution of kick 
velocities as a function of the symmetric mass parame- 
ter fitted nicely Fitchett's law [25] whose original deriva- 
tion was based on post-Newtonian analytical estimates of 
gravitational wave emission in inspiral binaries. This cor- 
responds to a Newtonian signature in the dynamics and 
this result may suggest that even in the head-on post- 
merger phase, described by RT dynamics, a component 
of the post-Newtonian dynamics of two interacting bodies 



emitting gravitational waves might be present/preserved. 
In the case of non head-on collisions examined in the 
present paper a similar Newtonian signature appears in- 
volving the relation between the incident angle of colli- 
sion Po and the scattering angle of the remnant /?/, as we 
discuss in Section X. 



IX. THE ANGULAR WAVE PATTERN AND 
THE BREMSSTRAHLUNG REGIME OF THE 
GRAVITATIONAL WAVES 

As discussed in Sec. II, the invariant characterization 
of a wave zone in RT spacetimes and the consequent 
presence of gravitational waves is given by the functions 
D{u,9,(l)) and B{u,9,(j)), Eqs. (|H])-(|ni), that characterize 
the curvature tensor components ([7]) at the wave zone. 
These functions correspond to the two modes of polariza- 
tion of the gravitational waves and contain all the infor- 
mation of the angular and time dependence of the gravi- 
tational wave amplitudes via the particular combination 



D(u,e,<f)) + iB{u,e,<f>) (rvE'4 



(72) 



where ^4 is the Weyl spinor associated with 
Nabcd/t 111, m. According to ©-(H]), the quantity 
([7^ is specified once we have the function A'(m, 9, (p) = 
1/P{u,9,(j)), which is in turn numerically obtained via 
the approximation (I35|) . 

In Fig. ini we display the polar plots of V D'^ + B"^ at 
early times u = 0.01 (dotted) ,m — 0.05 (dash-dotted) 
and u = 0.1 (continuous), and initial data parameters 
a = 0.2, Po = 55° and 7 = 0.5, with section by the plane 
(f) — Q° (corresponding to the plane of collision x — z). 
The plots in Fig. |9]show, for each time, a pattern with 
two dominant lobes in the forward direction of motion of 
the merged system, the first quadrant of the x — z plane. 
The direction of the Bondi momentum vector at m = 0.01 
makes an angle 85 ~ 8, 43° with the z axis. The pattern 
is typical of a bremsstrahlung process due to the decelera- 
tion of the merged system, analogous to the electromag- 
netic bremsstrahlung of a charge decelerated along its 
direction of motion. As time increases we observe that 
the cone enveloping the dominant lobes opens up and 
the amplitudes decrease. For later times (cf.Fig. [TT] for 
u — 5.0) the pattern evolves to the expected quadrupole 
structure with a much smaller amplitude. We mention 
that the increase of the initial boost parameter 7 would 
sharpen the forward cone enveloping of the two dominant 
lobes in the early regime, as expected in a ultrarelativis- 
tic configuration. In our computations we fixed mo = 10. 
In the Figures the z direction corresponds to the vertical 
axis, vertical axis). In Fig. [TU]we show the polar plot 
of y/lP+B^, with section by the plane (p = 90° (corre- 
sponding to the plane y — z, orthogonal to the plane of 
collision), at u ~ 0.1. The same initial data parameters 
of the previous Figure were used. As expected the pat- 
tern is symmetric about the z axis in accordance with 
the conservation of P^(u) = 0. We then see that, al- 
though gravitational waves are emitted outside the plane 
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FIG. 9: Polar plot of V-D^ + -B^ (section by the plane <?!) = 0°, 
corresponding to the plane of the collision) for times u = 0.01 
(dotted), u — 0.05 (dash-dotted) and u = 0.1 (continu- 
ous), and initial data parameters a = 0.2, po = 55° and 
7 = 0.5. The Figure shows a typical bremsstrahlung pattern, 
corresponding to a strong deceleration regime at early times, 
with two dominant lobes along the direction of motion of the 
merged system. The cone enveloping the two dominant lobes 
opens up as u increases. For larger values of 7 the envelop- 
ing cone of the forward lobes becomes more sharp. The z 
direction corresponds to the vertical axis. 
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FIG. 10: Polar plot of ^WTB^ (section by the plane 
(j) = 90°, corresponding to the plane y — z, orthogonal to 
the plane of collision) for a time u — 0.1, and initial data 
parameters a = 0.2, po — 55° and 7 = 0.5. The symme- 
try about the z-axis is in accordance with the conservation of 
P^{u) = 0. Although gravitational waves are emitted out- 
side the plane of collision, the zero net momentum flux of this 
radiation component is consistent with the planar nature of 
the collision. 



of the collision, this radiation component has a zero net 
momentum flux. Therefore it does not extract momen- 
tum of the system, consistent with the planar nature of 
the collision. 

Finally in Fig. [TT] we display the polar plot of 
y/D^ + B^ (section by the plane 4> = 90°) for a time 
w = 0.5 and for the same initial data parameters of Fig. 
O We can see the opening of the dominant lobes and the 
forming of the final quadrupole angular pattern. 



X. A NEWTONIAN SIGNATURE IN THE 
RELATION BETWEEN THE INCIDENCE 
ANGLE AND THE SCATTERING ANGLE OF 
THE REMNANT 

One of the important parameters of the remnant black 
hole obtained from our numerical experiments is the an- 
gle pf - defined a.s pf = cos~^(n3/) and denoted the 
scattering angle of the remnant - resulting from the evo- 
lution of initial data with initial incidence angle po . In 
general, for the mass ratio parameter a sufhciently small, 
Pf is much smaller than the incidence angle po. However 
for large values of a this behavior changes drastically as 
is the case for equal or nearly equal mass black holes. 

No simple analytical relation between po and p/ was 
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FIG. 11: Polar plot of ^fWT^ (section by the plane = 
0°) for a time u — 5.0, and same initial data parameters of 
the previous Figures, showing the opening of the lobes and 
the setting already of the final quadrupole angular pattern. 
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possible to be obtained from the Bondi-Sachs conserva- 
tion laws. However, as we will show, the scattering angle 
can be related to the incidence angle by a formula that 
closely approximates the Newtonian relation of angles in 
the nonelastic scattering of classical particles, and which 
depends basically only on the mass ratio a. To see this 
let us consider the schematic diagram of collision shown 
in Fig. [T21 There we have the two initial colliding black 
holes, one with total mass M boosted along the positive 
z axis and the other with total mass aM, with a < 1, 
bosted along a direction making an angle pq with the z 
axis, in accordance with the asymptotic result (j23p . In 
this picture we consider the approximation in which we 
neglect the contribution of the initial gravitational wave 
content to the individual masses, although this contri- 
bution is small and in any case can be isolated. The 
direction of the momentum of the remnant black hole 
makes an angle pf (the scattering angle) with the z axis. 
Furthermore, although in RT dynamics the presence a 
global apparent horizon is already present at the initial 
time and the system may be viewed as a deformed black 
hole, we are also assuming that inside the event horizon 
we still have two individual colliding black holes. 

Under the above assumptions and taking into account 
that the dynamics of the collision involves mass and mo- 
mentum loss due to the emission of gravitational waves, 
the balance of linear momentum can then be approxi- 
mately expressed as 

M sinh7 (1 — a cos po) = (Sz) mo A"/ cos p/ sinh7/, 
A4"sinh7 (asinpo) = {S^) moA'/ sin p/ sinh 7/, 

where Sz and dx are correction factors introduced to 
account for the mass and momentum loss in the process. 
The ratio of the above equations yield the relation 



AH 



a sm Po 
1 — Q cos po 



— S tan Pf 



(73) 



where 6 = (dx/Sz)- Eq. (|73|) provides a relation be- 
tween the initial incidence angle and the scattering angle 
of the remnant and must be validated by the results of 
the numerical experiments with the initial data (j53p . We 
note that d73l) reproduces the equation for the inelastic 
collision of classical particles when 6 = 1. 

Numerical results for pf obtained from initial data 
sampled in the whole range < a < 1 shows that the 
parameter i5 ~ 1, its discrepancy from the Newtonian 
classical value 1 occurring for relatively small values of 
a, where the asymmetry in the data is large. For illus- 
tration 

{a = 0.1, po = 45°, p/~4.16°, 5 ~ 1.046047) 
{a = 0.6, po = 45°, Pf ~ 35.78°, S ~ 1.022573) 
(a = 1.0, po = 45°, Pf ~ 67.50°, S ~ 1.000000) 

As a increases (5 — > 1, the limiting value 6=1 being 
attained for a = 1.0 and any po, yielding consequently 
(for a = 1) 




Pf = 



(180° - po) 



(74) 



FIG. 12: Schematic diagram of collision of two boosted black 
holes with initial infalling velocity parameter 7 and mass ratio 
a. The dashed line depicts the common apparent horizon 
enclosing the two colliding black holes. 



with error of the order of, or smaller than 10^^, for all 
experiments made with a = 1. In fact, pf = (180°— po)/2 
is an exact solution for the Newtonian inelastic collision 
equation (namely, Eq. ((75)) with 6 = 1). It is also worth 
remarking that for equal-mass initial colliding black holes 
{a = 1) the direction determined by the angle (180° — 
Po)/2 with the z axis corresponds to an axis of symmetry 
of the initial data (cf. Fig. [T2l) . 

Finally it can be of physical interest to compare the 
scattering angle with the initial momentum direction. In 
our above approximation of the momentum conservation 
we obtain, by a straightforward calculation, that the an- 
gle Oo of the initial momentum direction with the z-axis 
is given by tan9o = P'"(uo)/P^(«o) = ct sin po/(l— acos po), 
resulting from the relation ()73p that tanSo = S tan pf. 
This shows that the direction of motion of the remnant 
presents a very small deviation from the direction of the 
initial momentum as seen from an asymptotic iner- 
tial observer. This small difference is connected to the 
small values of the impulse imparted to the merged sys- 
tem by the gravitational waves emitted. However these 
small values of the impulse will be sufficient to generate 
kicks of a few hundred km/s in the system. 

The above results, relating po and pf, suggest a New- 
tonian signature present in the dynamics of the merged 
system with initial data (j53p , indicating that the individ- 
uality of the two initial colliding black holes is in some 
sense preserved in the initial deformed merged system. 



XI. CONCLUSIONS AND FINAL REMARKS 

In this paper we have examined the numerical evo- 
lution of characteristic initial data corresponding to non 
head-on collisions of two Schwarzschild black holes, in the 
realm of Robinson- Trautman non-axisymmetric space- 
times. The Robinson- Trautman spacetimes present al- 
ready a global apparent horizon so that the dynamics 
covers a regime where the merger has already set out. We 
have constructed initial data for the characteristic surface 
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formalism of Robinson- Trautman dynamics, that repre- 
sent instantaneously this system. The RT equation is 
integrated numerically using a Galerkin spectral method 
with a projection basis space in two variables (the spheri- 
cal harmonics). The planar nature of a generic non- head- 
on collision of two black holes, which is restricted to the 
plane determined by the two directions of motion of the 
two initial colliding black holes, is used to test the ac- 
curacy of our numerical method and allows us to save a 
lot of computational effort. The three basic independent 
parameters characterizing the initial data are the mass 
ratio parameter a, the boost parameter 7 (that defines 
the initial infalling velocity v = tanh7 of the initial in- 
dividual black holes) and the incidence collision angle po 
of the two initial colliding black holes. 

During the merger gravitational waves are emitted, ex- 
tracting mass and linear momentum of the system, until 
the remnant configuration is reached, when the gravita- 
tional wave emission ceases. Typically the remnant is a 
boosted Schwarzschild black hole defined by its final rest 
mass niQKj, its final velocity parameter Vf — tanh7y 
and its direction of motion determined by the remnant 
scattering angle pf. 

We evaluate the efficiency A of the mass-energy ex- 
traction as a function of the mass ratio parameter a for 
several values of the incidence angle po. We obtain that 
head-on collisions {po = 0°) constitute an upper bound, 
while orthogonal collisions (po = 90°) constitute a lower 
bound for the efficiency. This fact is connected to the 
behavior of the total gravitational wave energy emitted 
in the process since the head-on collision and the orthog- 
onal collision are configurations that correspond, respec- 
tively, to an upper bound and a lower bound for the total 
energy emitted. The efficiency presents a maximum at 
a = 1 for all values of the incidence angle, 0° < po < 90°, 
and for any boost parameter 7, which is shown to be a 
consequence of some symmetries of the initial data (I53p . 

The analysis of the momentum and energy flux of 
the gravitational waves is made using the Bondi-Sachs 
energy-momentum conservation laws. The momentum 
extraction by gravitational waves with the consequent re- 
coil of the system is evaluated, with respect to an asymp- 
totic Lorentz frame. For incidence angles po up to a cer- 
tain threshold value, we have that the net momentum 
flux carried out by the gravitational waves along the z 
axis is positive for a short initial phase leading to an in- 
crease of the momentum of the merged system along this 
direction, dP^{u)/du > 0. This phase is followed by a 
regime with negative net momentum flux that lasts un- 
til the final remnant configuration is attained. For po 
larger than the threshold value the initial phase of posi- 
tive momentum flux is no longer present. For a = 0.1 the 
threshold angle is po ~ 61°. Relative to the asymptotic 
Lorentz frame the Bondi-Sachs momentum components 
P^{u) and P^(u) are positive for all u. 

On the other hand, the net momentum flux along the a; 
direction is negative, for all times u and incidence angle 
poi so that the momentum of the system along the x 



direction always decreases. We have that P^^{u) = for 
all u, a consequence of the planar nature of the collision, 
so that the momentum of the merged system along this 
direction is conserved (in fact we have P^(u) = for all 
u). 

In general we have a net strong deceleration of the 
system due to the emission of gravitational waves, as a 
consequence of the the total impulse of the gravitational 
waves corresponding to the integrated fluxes. Typically, 
for large u ~ u/, the curves of the impulse tend to a 
negative constant value (a plateau) associated with the 
flnal configuration when the gravitational wave emission 
ceases. As a consequence, in accordance with the impulse 
conservation (j6ip for ^ — Uf, there is a net momentum 
decrease due to the total impulse of the gravitational 
waves imparted to the binary merged system. The im- 
portance of the initial acceleration of the z-component 
of the momentum (when present) is to produce a large 
inspiral branch in the motion of the center-of-mass of the 
system, as seen from the initial zero momentum frame, 
as we will discuss in a future publication. 

The momentum distribution of the remnant Pf/rriQ 
was evaluated for several po in the range 0° < po < 90. 
It has the head-on collision case as a lower bound and the 
orthogonal collision as an upper bound, for all a. A new 
feature is observed in the final momentum distribution 
for equal- mass black holes, in the case of non head-on 
collisions. Contrary to the case of a head-on collision 
and of the merging of black hole binary inspirals, for a 
non head-on collision of two equal mass black holes the 
net gravitational wave flux emitted is nonzero. There- 
fore the final momentum Pf for a = 1 is nonzero, for any 
Po ^ 0, even if measured in an inertial frame with velocity 
Vi„ = 'P{uq)/Mb{uq) (the zero initial-Bondi-momentum 
frame) relative to a rest inertial frame at infinity. The 
issue of kick velocities of the merged system due to the 
recoil of gravitational wave emission is currently being 
examined. Our analysis is based on the impulse conser- 
vation equations (j6ip and we obtain that, in a zero initial 
Bondi momentum frame, the velocity of the center-of- 
mass a,t u = Uf coincides approximately with the kick 
velocity. The distribution of the kick velocity as a func- 
tion of the symmetric mass parameter satisfies a modified 
Fitchett-Blanchet law[25j. This modification is necessary 
to give account of the fact that, in a non head-on collision 
of two equal-mass black holes, the net gravitational wave 
flux emitted is nonzero, contrary to the case of a head-on 
collision and of merging of black hole inspirals. 

The angular pattern of the gravitational waves emitted 
(with its two modes of polarization included) is examined 
for the case of a low mass ratio parameter and a large 
incidence angle. In the plane of the collision the pattern 
is typically bremsstrahlung, with two dominant lobes in 
the forward direction of motion of the merged system, 
characteristic of a strong initial decelerated regime of the 
merged system. We show that gravitational waves are 
also emitted outside the plane of the collision but that 
this radiation has a zero net momentum flux. Therefore it 
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does not extract momentum of the system, in accordance 
with the planar nature of the collision. 

Finally we also obtained a relation between the initial 
incidence angle po and the scattering angle pf that closely 
reproduces a result for the inelastic collision of classical 
particles in Newtonian dynamics, and is validated by the 
numerical results. For a = 1 we recover numerically the 
exact Newtonian result, with the p/ = (180° — po)/2 be- 
ing a consequence of the symmetry of the initial data as if 
the two initial black holes and the remnant were individ- 



ual classical particles. This result suggests a Newtonian 
signature present in the dynamics of the merged system 
with initial data ()53l 
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